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ABSTRACT 

The backward evolution approach to modeUing galaxy source counts is re- visited in 
the wake of the numerous results and revelations from the Infrared Space Observatory 
(ISO), the Submillimetre Common User Bolometer Array (SCUBA) and the detections 
and measurements of the cosmic extragalactic background light. Using the framework 
of the Pearson & Rowan-Robinson (1996) galaxy evolution model, the observed source 
counts and background measurements are used to constrain the evolution in the galaxy 
population. It is found that a strong evolution in both density and luminosity of the 
high luminosity tail of the IR luminosity function, interpreted as the ultraluminous 
galaxies discovered first by IRAS and later elevated in status by SCUBA and ISO, can 
account for the source counts from 15|xm (where it matches the undulations in the 
integral counts and hump in the differential counts extremely well), to the sub- mm 
region, as well as explain the peak in the cosmic infrared background at ~140|xm. 
The sub-mm counts are interpreted as the superposition of 2 separate populations 
comprising of ULIGs at the brighter sub-mm fluxes and starburst galaxies at fluxes 
fainter than ~2mJy. In this scenario the high rcdshift ULIGs are tenuously interpreted 
as the progenitors of today's giant elliptical galaxies. 

All the source count models can be accessed via the world wide web at this URL 



http : / /www. ir.isas.ac.jp/ ^ .pp/ counts / 



Key words: Cosmology: source counts - Infrared: source counts, background - Galax- 
ies: evolution. 



1 INTRODUCTION 

The formation epoch of galaxies has long been the holy grail 
of cosmology. Traditionally one of the most popular and suc- 
cessful ways of glimpsing this epoch has been the analysis of 
extragalactic source counts or more commonly galaxy num- 
ber counts. Although originally intended as a means of de- 
termining the geometry of the Universe this tool has shown 
to be much more viable and useful in the study of the evolu- 
tionary and st ar for mation history of galaxies (Kirshner et 
al.(1981), Ellis(1987)). In particular, the large dust content 



of galaxies discovered by IRAS and the strong K-corrections 
have made the infra-red and more recently, the sub-mm, ex- 
tremely important wavebands for the study of galaxy evolu- 
tion and cosmic star formation. By modelling the counts of 
galaxies it has become possible to analyze the evolution of 
the galaxy population out to redshifts of unity and higher. 



Further information contact Chris Pearson (cpp@ir.isas.ac.jp) 
http : / /www .iT.isas.ac.jp/ ^ cpp / counts / 



Coupled with the advent of new observational data, galaxy 
number counts are poised to discriminate between rivalling 
evolutionary theories. 

Recently, it has become the vogue to tackle the prob- 
lem of the number distribution of galaxies in the Universe 
via one of two metho ds. The first (of ten called the backward 
evolution approach (Lonsdale 199C)), takes the observed, 
present day [z = 0) luminosity function and evolves it in 
luminosity and/or density back out to higher redshifts as- 
suming s ome p arameterization of the e volution (Beichman 
& Helou ( |l99l|), Blain & Lo ngair ( p93| ), Pearson & Rowan- 



Robinson (L996), Xu et al. (1996), Takeuchi et al. (1998)) 



This method has the advantages that it is both direct and 
relatively simple to implement with few free parameters and 
assumptions about the Universe at earlier times. The disad- 
vantage in the past was that the information on which the as- 
sumptions about evolution were made often came from IRAS 
data which extended out only to relatively low redshifts 
~ 0.02 — 0.2. In the alternative method, known as the for- 
ward evolution approach the evolution is computed assuming 
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some initial conditions of the physical processes of chemi- 
cal evolution and photometri c evolution of the stella r pop- 
ulations that heat the dust ( Franceschini et al. 199^ ). The 
so-called semi analytical approach combines these assump- 
tions on the chemical/photometric evolution of galaxies with 
models considering the dissipative and non-dissipative pro- 
cesses affecti ng ga laxy formation within dark matter haloes 
(Cole et al. (1994), Baugh et al (1996)) and has provided a 



reason able fit to the source counts in the infra-red (Guider- 
doni et al. 1998). The forward evolution approach has the 



advantage of being based on perhaps a more fundamental 
set of assumptions, however the obvious disadvantages are 
the larger number of free/unknown parameters assumed in 
these models. 

Furthermore, with the recent adv ent o f large, deep 
survey s by SCUBA (Hollan d et al. ( [l999| ), e.g. Sma ll 
et al.( [l997| ) , Hug hes et al. (|l998|), B arger et al.( [l998| ), 




Blain e t al.(|l999a|) , Barger et al.(| l999)) an d ISO (Kessler 
et al.(1996), e.g. Bogun e t al. (1996), Serjea nt e t 
al. (yji97|). Tani guch i et al. ( |l997| ), Elb az et al. ( [l998| ), 
Stickel 
al. (pi99a[) , 

al. (|2000a| ), 

al. (2000b)), including constraints on the source counts 
to faint e r leve ls and the IR bac kgrou nd (Lagache et 



Kawara et al. 
Altieri et al. 
Lindern-V0rnle et al. 



(|1998D , Flores et 
fp,999|) , Oliver et 
(pOOOj), Oliver et 



al.(1999),( |2000[ ), Matsuhara et al. ( pOOOj ), Matsumoto et 
al.(2000)) and detections by the COBE FI RAS/ DIRBE in- 
struments at 140 & 240v im (P uget et al. (1996), Fixsen et 
al. (1998), Mauser et al. (1998)) the tools are now available 
to constrain the backward evolution methodology to signifi- 
cantly higher redshifts. In the light of these recent advances 
in the mm-IR region the backward evolution method is re- 
visited. 

The aim of this paper is not to produce a definitive 



evolutionary model to fit all observed counts while satisfy- 



ing various preferences on the various dominant populations 
(starburst, AGN, luminous infrared galaxies (LIG's Ljr > 
IO^^Lq), ultraluminous galaxies (ULIG's L/h > IO^^Lq), 
hyperluminous galaxies {Lm > IO^'^Lq)). Rather I seek to 
investigate the various avenues of evolution that can viably 
fit the galaxy source counts from the mid-IR to sub-mm 
wavelengths whilst not violating the constraints set by the 
IR background and CMB measurements. The paper is struc- 
tured as follows, in section ^ the galaxy count and world 
model parameters are outlined. Section ^ reviews current 
evolutionary models and constraints and analyszes simple 
density and luminosity evolution scenarios. Section ^ for- 
mulates a new evolutionary source count model and applies 
it to the observed source counts and background measure- 
ments from the mid-IR to sub-mm wavelengths. The possible 
interpretations and conclusions are reported in sections ^ & 
^ respectively. 



2 MODEL PARAMETERS 

In general, for observations at frequency v = cjX, the total 
number of sources per steradian observable down to a flux 
sensitivity S is given by; 



iV(S.) = , 



(1) 



dV{z) 
dz 



fi-o(l + z)2(l-ha„«)l/2' 



(2) 



where — d^/dlgL — the zero redshift differential lu- 
minosity function per decade in luminosity L, parameteriz- 
ing the number density of extragalactic sources as a func- 
tion of luminosity at a frequency z/. dV/dz is the general 
form of the differential volume element with redshift, z, and 
is dependent on the luminosity distance, D, and assumed 
cosmological world model. f{z) & g{z) are luminosity and 
density evolutionary parameters respectively. The limiting 
redshift z{L,S), is determined by the flux, S, below which 
a source of luminosity, L, becomes too weak to be included 
in the sample of galaxies observed, where; 



v.L 



47r_D2 UoL 



D{z) 



2c 



{n,z + [Qo - 2)[(i + nozf" - 1]), 



(3) 



(4) 



where 



^ is effectively the ratio of the emission lu- 
minosity (luminosity at Pe) to the observed luminosity (i.e. 
the K-correction). Even assuming a single luminosity func- 
tion and evolutionary history for all extragalactic sources 
can in fact come impressively close to matching the true 
numbers and distribution of sources over a wide range of 
wav eleng ths as o bserved in the Universe (e.g. Takeuchi et 
al. (|1999[), (^OOOp). 



Pearson & Rowan- Robinson (199(:) (hereafter cppRR) 
provided an impressively solid framework from which was 
modelled the source counts from radio frequencies, sub-mm, 
IR, optical, UV wavelengths to X-ray energies. The origi- 
nal IR models of cppRR consisted of a 4 component pa- 
rameterization of the extragalactic population based on the 



IR co lours of IRAS galaxies (Rowan- Robinson & Crawford 



198£). The model incorporated a normal galaxy (IR cirrus 



dominated) component defined by cool 100/60^Lm colours, 
a starburst component defined by warm 100/60 (i.m colours 
extending to an ultra/hyperluminous component at higher 
IR luminosities and an AGN Seyfert/QSO 3-30ii.m dust 
torus component defined by low 60/25nm colours. The basic 
framework of this model has been maintained in this further 
analysis. 

The original template spect ral energy distributions 
(SED's) of Row an- Robinson ( |l992| ). Rowan- Robinson & Ef- 
stathiou (199J) for the coo/ normal, uiarm starburst & ultra- 
luminous galaxies are replaced by the new radiative trans- 
fer models of Efstathiou & Siebenmorgen (2000) and Ef- 
stathiou. Rowan- Robinson & Siebenmorgen (2000) respec- 
tively. The new starburst galaxy models are calculated con- 
sidering the evolution of an ensemble of optically thick giant 
molecular clouds, illuminated by embedded massive stars. 
The starburst models are defined by the parameters t , the 
age of the starburst in Myrs; t, the initial optical depth (in 
y ) of the molecular clouds; Xi the ratio of radiation field 
to the local solar neighborhood. The starburst SED can be 
considered akin to M82, while the ULIG SED represents 
an ARP220 type galaxy. The introduction of this ARP220 
type SED replaces the more extreme h yperlu mi nous galaxv 
IRAS F10214 (Rowan-Robinson et al. ( |l991a|) , ( |l993| )) that 
was used in the earlier models of cppRR. F10214 does not 
have a typical ULIG SED and such a hyperluminous galaxy 
population may represent an even more extreme extension of 
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the starburst p henomenon than the U LIGs, only 39 being so 
far identified ( Rowan- Robinson 2000 ). Efstathiou, Rowan- 
Robinson & Siebenmoreen ( |200q ) demonstrated that ex- 
ponentially decaying 10^ ~ lO^yrs old starbursts are capa- 
ble of reproducing both the IRAS colours of galaxies and 
also provide good fits to the data for the stereotypical star- 
burst galaxy M82 a nd to the recently I SO observed starburst 
galaxy NGC6090 (lAcosta et al. 19961). 



Of particular importance to the models in the mid- 
infrared region is the inclusion of the unidentified infra-red 
bands (UIB), popularly supposed to be band emission from 
Pol ycyc lic Aromatic Hydrocarbons (PAH's (Puget & Leger 



198!;)), into the model SED's. T hese features ar e prominent 
from 3-12|jjn i n most galaxies ( (Lu et al. 1997 ), (Boulade 
1999), (Vigroux L. et al. 1996D) having equivalent 



et al. 

widths of as much as 10 (xm, comparable to those of the ISO- 
CAM MIR filter band passes (3.5^m & 6^m for the LW2(5- 
8.5ii.m) filter and LW3(12-18(i.m) filters respectively). These 
UIB's are expected to significantly distort the mid-IR source 
counts compared to those cons tructe d purely on the basis of 
the IRAS SED data (Xu et al. (|l998[), Pearson et al. (|200l|)). 



A further significant and necessary improvement over the 
original cppRR model is the introduction of t he full range of 
PAH features ( Biebenmorgen fc Krugel 1992 ) that were not 
considered in the earlier model predictions. Fig. |l] shows the 
new model SED's. The shape and form of the PAH features 
assumed ii i diffe rent models can v ary significantly (e.g. see 



Xu et al. (1998), Dole et al.(200C)) but although the PAH 



features may seem overly strong or sharply defined, as long 
as the energy contained within the feature peaks is consis- 
tent with observation and the SEDs are correctly smoothed 
and convolved with the filter response function, the general 
profile of the individual features is not critical as galaxies 
will naturally have a range of feature profi les in their SEDs. 

Note that recently Dunne et al. ( 200C ) have shown from 
SCUBA data on local galaxies that there is a correlation 
between the far-IR-sub-mm SED and the 60^.m luminosity 
function. This effect is not treated in the model SEDs and 
will be the subject of a later investigation. 

For the AGN SED the original SED of Rowan- Robinson 



( [1995| ) is retained. The MIR SED of AGN is relatively fea- 
tureless compared to that of starburst galaxies, exhibiting 
none of th e broadband UI B features present in the latter's 
spectrum (Lutz et al. 1997). The AGN central engine (black 



hole) would be expected to produce an extreme environ- 
ment destroying or changing the smaller dust grains and 
their composition. 

The model components are distinguished on a basis of 
their SED's and number-luminosity distributions (i.e. lumi- 
nosity function) . The cool normal galaxies, warm starbursts 
& ULIG galaxies are well described locally by the IRAS 
60jxm luminosity function (saunders et al. 1990) which is 



a log Gaussian function producing a slightly broader distri- 
bution than a traditional Schecter function and preferable 
to the slightly l ess physical powe r law parameterizations of 
IRAS galaxies (Soifer et al. 1987). In this parameterization 



the normal and starburst galaxies extend in 60 (xm luminos- 
ity from IgLaofim = 8 — 11.5Lq while the ULIG population 
forms the high end tail of the starburst luminosity function 
from IgLeoum ~ 11.5 — ULq. The AGN are represented 
by the 12 ^.m luminosity function of Lawrence et al. (198(: ) 
defined from lgL\2ij.m = 8 — 14^0 with parameters cor- 




12.5 



13 13.5 
Ig (V) {Hz} 



14.5 



Figure 1. Model spectral energy distributions used for the source 
count models. 



responding to the 12)i.m sample of Rush et al. ( 1993| ) for 
Seyfert 1 and Seyfert 2 galaxies. The parameters for the dif- 
ferent components given in table ^ correspond to those in 
the equations below for the 60nm galaxy and 12|i.m AGN 
luminosity functions respectively; 



0(^6 



0(Li 



d£ 
dZ 



exp[ 



1 

'2a^ 



dL 



'(1 + 



L*l3 



^13 



(5) 



(6) 



To shift the luminosity function from the wavelength 
at which the luminosity function is defined, Xlf, to the 
observation wavelength, Xots, the ratio L(\obs) / L{Xlf) is 
obtained via the model template spectra. 

The model assumes n = 0.1, Ho = 50fcms-^Mpc-\ 
Alternative cosmologies are not examined in this work so 
as to ensure that a clear distinction between the evolution 
scenarios can be made. The maximum redshift assumed for 
the models is z = 7 although assuming a higher redshift of 
10 makes little difference to the results presented here. 



3 EVOLUTION 

3.1 Observational Evolutionary Evidence 

The original cppRR model assumed a simple power law lu- 
minosity evolution L{z) = L{z = 0)(1 -I- z)^'^ for the star- 
burst, ultraluminous and AGN co mponents and was found 
to fit the IRAS data extremely well (lOliver et al. 19921). This 



evolution was motivated by similar evolutionary models at 
radi o and X-ray wavelengths (Boyle et al. (1988), Benn et 
al. (1993), Condon (1994)). Recently, similar evolution has 



also been obse rved in optically selected starburst galaxies 
( Lilly al. 1996 ) and in the submilliter emissio n from radio 
loud gal axies observed at 850[J.m by SCUBA (Archibald et 
al. 2000). Although this simple form of evolution still pro- 



vides a satisfactory fit to source counts down to quite faint 
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Table 1. Luminosity function parameters. 



Component 


a 


a 


/3 




C 


cool 60 [am 


1.15 


0.463 




9.62 


2.000 X 10-3 


warm 60|-im 


1.27 


0.626 




9.99 


3.250 X 10""' 


ULIG 60 (om 




0.1 




11.6 


4.185 X 10-^ 


Seyfert 1 12 ^m 


1 




2.1 


9.552 


2.0969 X lO""* 


Seyfert 2 12^m 


1 




2.5 


9.952 


8.0251 X 10"^ 



sensitivities (ImJy at IS^tm, Serjeant et al. ( ^000| ) , ~0.3Jy 
& ~0.1Jy at 170ii-m & QOjtm respectively), it cannot match 
the observed bumps and dips in the ISO 15vim counts below 
ImJy and falls below the observed faint counts at far-IR 
wavelengths (e.g. 170n.m FIRBACK; P uget e t al. short- 
citepug99, Lockman Ifole; Kawara et al. (1998)). Moreover, 
the extremely high counts reported by the SCUBA instru- 
ment at submillimetre wavelengths at fluxes brighter than 
~2mJy cannot be reproduced in any way by the cppRR 
luminosity evolution model, although it should be noted 
that the predictions are tantalizingly close to the observed 
counts at fluxes fainter than Sssoum < 2mJy. Furthermore, 
although the faintest source counts from the60p.m IRAS 
sur vey (V FSS, ~120mJy Gregorich et al. ( |l995| ), Oliver et 
al. (1995), Berlin et al. (1997)) could be well explained by 
the cp pRR pure evolution model at these relatively bright 



cally (2 < 0.1), ULIGs have a space density approximately 
similar to optically selected QSOs, the analysis of the IRAS 
FSC hinted at much stronger evolution at higher redshifts to 
z « 1 thus outnumbering the QSOs by a factor of 2. Based 
on the premise that the ULIG pop ulation may follo w a sim- 
ilar evolutionary track as QSO's (Hassinger 1998), Kim & 
Saunders ( 



1995 



modelled the ULIG population with den- 
sity evolution of the form (l-\-zy , with g ~ 7.6±3.2. Similar 
strong density evolution (equivalent to g ~ 10 at 2 ~ 1) has 
also been propose d to fit the deep 15|J.m ISO source counts 
(Dole et al. 200C). The vast majority of ULIGs appear to 
be in the la tter s tages of merging (98% in the sample of 



Borne et al. (199S;), see also, A rribas et al. ( ^000| ), Borne et 



al. ( |2000[ ), Colina et al. ( |2000| )) with double nuclei, overlap- 
ping disks and tidal fe atures clearly present in m a ny images 
(ganders et al. 199^ ) ( [Melnick fc Mirabel 1990| ), ( [Clement ' 



fluxes 1 he evolutionary constraints are extremely tenuous et al. 199(^ ) . The large amounts of dust present in the central 



with many different scenarios being allowed ( the IRAS sam- 
ples not probing to si gnificantly large redshifts (QDO T sam- 
ple Zmedian ~ 0.03 ( Rowau-Robiuson et al. 1991b), VFSS 
analysis z < 0.3). At redshifts higher than ~ 0.25, the evo- 
lution is not well constrained. 

Recently, several authors have attempted to combat 
these s hortcomings bv incorporating stronger or more dras- 



tic evo ution scenarios into the source count framework. 



Rowan- Robinson (200C) included evolution explicitly into 

the TT? firinm liiminnsity fTinrtron rp-prriHTiring many of 



the ob aervables if a low densitv cosmology was assumed 
(fio ~ 0.3). For higher values of the density parameter it was 
concluded that a new, strongly evolving popul ation would 
have t o be i nvoked (e.g. Franceschini et al.(1998)). Takeuchi 
et al. (200C) treated a global luminosity evolution scenario 
in a stepwise, non-parametric form concluding that strong 
evolution to z~l was necesscary to fit the ISO source counts 
at 15iim and 170^im along with the far-IR background. An 
advantage of the formulation of cppRR was that the 60^im 
luminosity function was treated as a superposition of galaxy 
populations providing the added fiexibility of specific evolu- 
tion with extragalactic class. However, these evolution sce- 
narios must be reviewed if we are to hope to fit the SCUBA 
counts at 850iim and ISO 15iim & 170p.m faint counts. 

From a statistical analysis of the 60^im IRAS VFSS, 
Berlin et al. (1997) concluded that the most luminous 
IRAS galaxies should be ~ 5'^7 times more numerous at 
z«0.3 than at the present epoch with IR luminosities of 
the order of Leo ~ 2x10^^ and optical magnitudes between 
20 < Bj < 21 with red colours {Bj — Rf ~ 1.6 making 



them potential U LIG c andidates (Clements et al. 199f))). 
Kim & Saunders (1998 ) compiled a complete fiux limited 



sample of ULIGs from the IRAS Faint Source Catalogue 
(FSC) with a flux /go > IJy to a maximum redshift of 
z — 0.268 and mean redshift < z >~ 0.15. Although lo- 



rcgions of these galaxies provides an ideal environment for 
circumnuclear starburst and/or AGN adding to the grow- 
ing evidence that ULIGs could be the progenitors of today's 
giant elliptical galaxies formed by merger induced collapse 
of the galaxy core, formation of globular clusters and galac- 
tic superwinds resulting in a dust free core an d correspond- 
ing enrichment of the intergalactic medium (Kormendy & 



Sanders 1992), (Surace et al. 199J). However, in the local 



Universe z <~ 0.08, the space density of ULIG s is approxi- 
mately .001 per sq.deg. or « 1.1 x Mpc~^ ( Saunders et 
al. 1990), a factor of as rnuch as «1500 lower than early type 



galaxies (Trentham 200C). Thus if ULIGs are to be present 
in sufficient numbers to be representative today's gE galaxy 
population then their space density must increase drastically 
with look back time (a factor of 1500 lO'^'Afpc-^). 



3.2 Evolution of the high luminosity tail of the IR 
luminosity function 

Under this motivation, extreme number evolution of the 
ULIG population is considered, i.e. the high luminosity tail 
of the 60nm luminosity function. Fig.^ shows the effect of 
such evolution on the 60p.m differential luminosity function. 
Evolution in the number of galaxies with redshift (density 
evolution) produces a vertical shift in the luminosity func- 
tion plane. Evolution in the luminosity of galaxies (luminos- 
ity evolution) produces a horizontal shift in the luminosity 
function plane. In reality a combination of the two scenar- 
ios may well be realized. In the figure, only evolution in 
the high luminosity tail of the luminosity function is shown. 
Here I have defined a population of luminous irffrared galax- 
ies (LIG/ULIG) with IgL* ~ 11.6Lq, i.e. corresponding to 

11.8Lq such that the 



a far-IR luminosity (50-100[J.m) of - 
total IR luminosity (l-lO OOum) of • 



12L0 (Hclou, Soifer 



fc Rowan- Robinson 1984). The definition of L* is given in 
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Figure 2. Evolution of the high luminosity tail of the 60 ^un dilfer- 
ential luminosity function. The zero redshift (no evolution) 60|jm 
luminosity function and high luminosity tail are shown in com- 
parison with an evolved L* = IO^^'^Lq luminous infrared galaxy 
luminosity function tail at z=l, shown for mild pure luminosity 
evolution, equivalent to (1 -|- z)*^'^, strong pure density evolution, 
equivalent to ~ (1 -|- z)^ and a combination of both. Luminosity 
evolution corresponds to a horizontal displacement in luminos- 
ity, density evolution corresponds to a vertical displacement in 
number density. 



equations 5 & 6 which define the luminosity functions used 
in these models. Although the luminosity functions are not 
classical Schecter functions, the definition of L* is almost 
identical. Changing the value of L* to higher luminosities 
shifts the luminous infra-red population down and to the 
right (towards a lower normalization) in the luminosity func- 
tion plane and has a similar effect on the corresponding cal- 
culated source counts as shown in figure ^ where a brighter 
L* results in detections at brighter sensitivity limits but 
at a lower number density (i.e. more powerful but scarcer 
sources). It should be kept in mind that the main aim of this 
work is to try to highlight the constraints, by means of sim- 
ple evolutionary scenarios, on the evolution of the high end 
tail of the 60(i.m IR luminosity function (luminosity func- 
tions in other wavebands are also avalible e.g. recently the 
sub- m m lum inosity function has been calculated by Dunne 
et al. (|2000[) using the SCUBA Local Universe Galaxy Sur- 



vey, SLUGS). 

As a starting point, I compare the pre-ISO view of 
galaxy evolution with the galaxy source counts observed at 
15, 60, 170 & 850nm by ISO, IRAS, ISO and SCUBA re- 
spectively. Three scenarios are presented for the evolution of 
the ULIG component, (a) no-evolution, (b) luminosity evo- 
lution of the form f{z) — {I + z)"^'^ (cppRR model) and (c) 
densit y evol ution of the form g{z) = (1 -I- 2)^ (e.g. Oliver 
et al. (1995)). The evolution is continued out to a redshift 



of 2.5 where it is assumed to level off and eventually de- 
cline t owards higher redshift. Fig. ^ shows the individually 
evolvin g ULIG components plus the"total source counts cor- 
responding to each scenario. The total source counts include 



the ULIG component and the normal, starburst, AGN com- 
ponents (evolving as L{z) = f{z)L{0) as per the cppRR 
model framework). From this simple representation of the 
evolution we can already ascertain important boundary con- 
ditions on the evolution of the ULIG component. At the 
brightest fluxes all models agree on the same normaliza- 
tion as they should seeing that at these fluxes the counts do 
not probe deep enough to discriminate between the different 
evolution scenarios. Both the luminosity and density evolu- 
tion scenarios are monotonically increasing in power with 
redshift so we may interperet a deficit with respect to the 
observed counts as a lack of evolutionary amplitude at lower 
redshifts while an excess with respect to the observed counts 
may be interpreted as an indication of too much evolution 
at higher redshifts. 

At 15(i.m, the integral counts appear to be fitted well by 
all scenarios with the only discrepancy between the density 
and luminosity evolution models appearing at the faintest 
fiux levels (remembering that evolution is also present in the 
other more est ablish ed components of the cppRR model). 
Serjeant et al. (2000) found that the 15(i.m integral counts 



could indeed be well fitted by a v ariety of evolving model s 
(cppRR, Franc eschini et al.(1994), Guideroni et al. (1997), 
Xu et al.(199^)). However, examining the 15iJ.m differen- 



tial counts provides a more discriminating picture. The dif- 
ferential counts effectively measure the slope of the inte- 
gral counts and are hence more sensitive to subtle changes 
that may not be prominent in the integral counts at first 
sight. Here, again we see a discrepancy between the two 
evolution scenarios with the density evolution model pro- 
ducing slightly higher differential counts (i.e. steeper inte- 
gral counts) at the faintest fiuxes (perhaps hinting at being 
marginally inconsistent with the observed counts due to too 
much power at higher redshifts?). However, what is more 
striking is that neither the classical density nor luminos- 
ity evolution models can reproduce the subtle changes of 
slope in the 15jim counts that are not easily apparent from 
the integral source counts alone. It is apparent from the 
differential counts that the counts steepen drastically from 
~3mJy-0.4mJy {lg{dN/dS) = algS,a « -3) above Eu- 
clidean expectations, turn over at ~0.4mJy and rapidl y flat - 
ten at fainter fluxes. The analysis by Serjeant et al. (200C) 
concluded that models incorporating pure luminosity evolu- 
tion provided a better fit than the density evolving models 
although they assumed a lower value for the density evolu- 
tion and moreover their observations did not extend down 
in detail below lOmJy in the differential counts. 

At 60ii-m the IRAS data does not really extend to deep 
enough fluxes to gain meaningful insight from the models, 
with both the luminosity and density evolution scenarios 
fltting both the integral and differential counts, although the 
faintest differential counts would imply stronger evolution 
than currently assumed. 

However, at 170nm none of the models can account for 
the observed source counts at fainter fluxes (although the 
brightest points are fltted by all models). At 170^.m we could 
envision a general increase in the magnitude of the evolution 
as generally fltting the observed counts (i.e. increasing the 
density evolution power law index g). However, in th e case 
of density evolution, fltting the FIRBACK counts ( Puget 



et al 1999) would result in violating the limits set by the 



fluctuation analysis in the Lockman Hole of Matsuhara et 
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Figure 3. Integral source counts of luminous/ultraluminous IR galaxies at (a) 15|^m, differential counts at (b) 15\A.m and integral counts 
at (c) 170|am, assuming values of IgL* = 11.2, 11.5, 12Lq. The values of Cquoted are the corresponding normalizations of the luminosity 
function given in equation |5| and table |^ . The shift with increasing IgL* to the lower right in the IgN/lgS plane reflects the higher 
brightness but lower number density (normalization) of the sources. 



al. (2000) if tlie trend of a flattening in the source counts is 
continued belo-w fluxes of 70mJy. This again points to some 
constraint on the magnitude of the allowed density evolu- 
tion if the evolution is allowed to increase monotonically to 
higher redshift. If density evolution is to account for the 
170vim counts while simultaneously not violating the fluctu- 
ation analysis limits then it must take place a relatively low 
redshif t (z~l). 



Fi nallv. in the sub-mm at 850fi.m. the biggest deficit be- 



tween the observed and predicted counts is se en. Surveys of 



bot h blank fields and around lensed clu sters ( Hughes et al 
199^)71 (^mafl, Ivison fc Blain 1997| ), ( |Barger et al. 1998} ) 



(Eales et al. 1999) of areas ranging from 0. 002-0. 12sq.deg., 
from l-8mjy have revealed source densities up to 3 orders 
of magnitudes above no-evolution predictions. Furthermore 
due to the steep K-corrections in the sub-mm, the fiux of 
high redshift galaxies is enhanced, resulting in almost no 
difference in the ability to detect a galaxy between redshift 
1-10. However, even at 850 (xm we also see that monotonically 
increasing power law density evolution cannot fit the counts 
and prod uces too much evolution at higher redshift. Blain et 
al. (1998) concurred with this result by finding that fitting 
the counts with pure density evolution would result in a vio- 
lation of the background radiation by 2 orders of magnitude. 
Interestingly however, is the luminosity evolution model that 
comes tentatively close to the observed source counts at the 
brightest and faintest fiuxes. 

From fig. ^ it is apparent that significantly drastic evo- 
lution, in both magnitude and form, will be required to fit 
the observed source counts. From the fits to the faint ends 
of the 170^im, 850nm and marginally to the ISjim counts, 
significant density evolution at high redshifts seems unlikely. 
As well as exceeding the observed source counts at 170ii-m 
& 850vim, such density evolution would also violate the con- 
straints set by the FIR background radiation. Matsuhara 
et al. (200C) found that the integral counts at 170ti.m in 
the Lockman Hole at a flux of ~150mJy had a steep slope, 
lg{dN/dS) — algS,a ~ —3. However they summized that 
if this slope did not become shallower before ~lmjy then 
the power in the calculated fluctuation spectrum would ex- 



ceed the observed power by a factor of 10. Therefore at some 
point the source counts must flatten to a value of q < —2. 
Fig. ^ shows that assuming monotonically increasing density 
evolution to z~ 2.5, the counts remain steep down to fiuxes 
fainter than ImJy. At 850 ^im, constraints set by COBE mea- 
surements of the sub- mm ba ckground radiation (Fixsen et 
al. ( |l998|) , Hauser et al.( [l998|) ) require that the SCU BA sub- 
mm count s i nust converge at fiuxes ~ 0.3 — 0.5m Jy ( Hughes 



et al. 1995), (Hughes 1999). Again, assuming power law den- 
sity evolution, the counts remain steep down to fiuxes deeper 
than O.lmJy. Observations and follow ups of the ISO 15jim 
surveys of the CFRS, HDF, Lockman Hole & Marano fields 
identified the sources with galaxies at redshifts of betwee n 
0.3 and 1 with a media n redshift ~0.8 (Flores et al. 1999a), 



( |Mann et al. 1997D , ( |Elbaz et al. 1999| ). 

In summary, there is too little density evolution at low 
redshift {z ~ 1), and too much density evolution at high 
redshift (2 ~ 2) in the power law evolution parameteriza- 
tion to fit the observed source counts at far-IR and sub-mm 
wavelengths. Furthermore, the 15ii-m and 170(i.m counts re- 
quire extreme evolution from a redshift of ~0.2 (the limit of 
IRAS) to ^l. On the other hand the SCUBA 850n,m counts 
and the 170ii-m fiuctuation analysis will constrain the up- 
per redshift limit to the evolution. Therefore, I assume that 
the ISOCAM 15(i.m differential counts can provide a good 
boundary to constrain the evolution in the low redshift Uni- 
verse (0 < 2 < 1) and the 850ii-m SCUBA sub-mm counts 
will constrain the evolution in the high redshift Universe 
(1 < 2 ~ 3). The evolution of the galaxy population will 
also be constrained at high redshifts by the measurements of 
the cosmic infrared background (CIRB). By simultaneously 
examining evolution scenarios that successfully fit both the 
850ii-m and 15jim counts it is expected that the model fits 
to the data at intermediate wavelengths and CIRB will also 
converge. 



3.3 Density Evolution Analysis 

Treating only the ULIG component, to investigate the den- 
sity evolution I consider 3 simple Gaussian models (see 
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Figure 4. Comparison of source counts including the luminous/ultraluminous IR galaxy (ULIG) component assuming the cppRR 
framework for the normal, starburst and AGN components. Integral and differential counts at 15|jm, 60|^m, 170|-mi & 850|-mi are shown 
for 3 classic evolution scenarios of the luminous/ultraluminous IR galaxy component {thin lines). Also shown is the total contribution of 
all components including the ULIG component (thick lines). A non-evolving ULIG component {solid lines). Pure luminosity evolution 
of the form f{z) = (1 + z)^'^ (Levol - dashed lines ) and pure density evolution of the form g{z) = (1 + z)^ {Devol - dotted lines). Data is 
from 15^m - Oliver et al. (1997), Aussel et al. (1999), Serjeant et al. (2000), Elbaz et al. (1999). 60nm - Lonsdale et al. (1990), Hacking 
& Houck (1987), Rowan-Robinson et al. (1990), Saunders (1990), Gregorich et al. (1995). 170|am - Kawara et al. (1998), Puget et al. 
(1999), Matsuhara et al. (2000) (fluctuation analysis - dotted box) . SbOytm - Small et al. (1997), Hughes et al. (1998), Bales et al. (1999), 
Barger et al. (1999), Blain et al. (1999). 
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Fig. ^, although not exponential in the strictest sense, here- 
after referred to as exponential evolution). Narrow Gaus- 
sian profiles peaking at a redshift of 1 and 2 are compared 
with each other and also with a scenario, referred to as 
the G aussian plateau, similar to the anvil model of Blain 
et al. (1999b), that evolves exponentially to redshift ~1 and 
declines exponentially after a redshift ~2. The Gaussian pro- 
files have the form, 



D( 



■ gexpi 



2cr2 



(7) 



with the peak magnitude of the evolution, g = 1000, the 
Gaussian width (effectively the evolution timescale) a = 0.1 
and the peak redshift = 1 & 2 respectively. The Gaussian 
plateau model assumes an exponential increase and decrease 
as equation ^, to and from redshifts of 1 and 2 respectively 
and a constant value of 1 + <? in between. Although these 
models are merely intended to highlight the effect of density 
evolution at different redshifts, the narrow Gaussian pro- 
files may be considered as epochs of short, intense merging 
while the Gaussian plateau represents a longer episode of 
galaxy merging (and consequent star formation / galaxy IR 
emission) . 

Using these density evolution scenarios, the correspond- 
ing evolving ULIG components are re-plotted in fig. ^ at 
15p.m and 850^im. Several points are worth noting in com- 
parison with fig. 1^. Assuming the exponential evolution from 
equation ^, a value of g = 1000 at redshift 1 corresponds 
roughly to a 3 ~ 10 in the (1 -I- z)^ power law scenario 
with the significant difference that the high redshift density 
evolution has vanished. The ISjim counts remain relatively 
unchanged since they are not so sensitive to any density evo- 
lution at high redshift {z > 1). However, the effect on the 
850 H-m sub-mm counts is significant in that the excess at the 
faint fluxes (~ ImJy) has disappeared. In the case of the 
Zp = 2 scenario, the ISjim counts cannot be fitted due to 
the fact that there is no evolution at lower redshift where it 
is needed. In fact, as can be seen from the Gaussian plateau 
model, the IS^im counts in general are clearly rather insen- 
sitive to density evolution at higher redshifts {z > 1). From 
the 15n.m counts, it is clear that power in the local density 
evolution is required. For the 850(i.m counts (fig. ^), density 
evolution peaking at either z=l or 2 is acceptable. However, 
the Gaussian plateau model, with constant density evolution 
from z=l-2 violates the faint end of the 850|i.m counts. 



3.4 Luminosity Evolution Analysis 

Despite the apparent success at 15vim, the brighter SCUBA 
counts at SSOjim (5 > ImJy) seem impossible to fit with 
this form of density evolution alone. Even the evolved counts 
still fall too far to the left and low in the lg{N) /lg{S) plane 
when compared to the observations > ImJy. In general, 
the brighter SCUBA population is comprised of fewer but 
more powerful sources (a high redshift 8mJy SCUBA source 
being approximately 10 more luminous than the most lu- 
minous local sources). Additional clues come from a further 
inspection of the luminosity evolution models of fig. ^. In 
general the addition of luminosity evolution has the generic 
effect of shifting the source counts towards brighter fluxes (to 
the right of the lg{N) /lg{S) plane (as the brighter, evolved 
sources are seen at lower sensitivities i.e. higher flux values). 



4 7 
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Figure 5. Density evolution parameterizations for the ULIG 
component models discussed in the text. Density evolution is in- 
vestigated via three exponential models. Two Gaussian parame- 
terizations peaking at a redshift of 1 and 2 respectively plus an 
exponential plateau model ranging from z=l-2. Also shown for 
comparison are classical power law density evolution models with 
power indices g=10 and g=7 (i.e. that previously assumed for the 
classical density evolution of IRAS sources, see fig. 4) respectively. 
The magnitude of the density evolution in the zp=l Gaussian sce- 
nario is equivalent to the magnitude of the evolution in the power 
law model using a value of g=10 at z=l. 



Therefore, I extend the analysis of the density evolution 
scenarios to include a similar subset of luminosity evolution 
models. Following the same doctrine as the density evolution 
analysis I define the luminosity evolution as. 



L{z) = 1-1- kexp[ 



2cr2 



(8) 



with the magnitude of the evolution, k = 5, a = 0.3, 
Zp — 1 &i 2 respectively. Thus this form would be equiva- 
lent to a value of fe ~ 2.3 for the case of {I + z)'' power 
law luminosity evolution. Note the disparity in the magni- 
tude / width of the luminosity and density distributions. 
It should be noted that, firstly, the luminosity and density 
evolution enter the cosmological equations differently. Den- 
sity evolution increases the number of sources observed in 
a given volume defined by a sensitivity S{z) and redshift, 
z. The volume remains the same and the number density 
of sources increase. On the other hand the effect of lumi- 
nosity evolution is to increase the observable volume via a 
lowering of S{z) while preserving the source number density. 
Furthermore, the luminosity evolution analysis is performed 
on the assumption of there already being the presence of 
density evolution at low redshift manifested as the Zp = 1 
exponential density evolution scenario. 

Fig. 1^ shows the differential and integral counts of the 
ULIG component at 15n.m and 850n-m respectively for the 
3 luminosity evolution scenarios described by eqn. pi In this 
analysis, density evolution of the form in eqn. M is assumed 
with parameter values of <; = 500, a = 0.1, Zp = 1. As 
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redshifts higher than unity allowing the Zp = 2 scenario to 
compete with the combined effect of the luminosity and den- 
sity evolution in the Zp — 1 scenario. However, it should be 
noted that there is a significant difference between the two 
Gaussian scenarios and the plateau scenario, ft would seem 
that constant luminosity evolution over the entire redshift 
range z=f-2 is ruled out by the count constraints. 

As with the density evolution analysis, the SCUBA 
850n-m counts are constrained by the evolution at higher 
redshift while the 15^,m counts are severely constrained at 
lower redshift. This bipolar constraint is highlighted in the 
way that the plateau curve hugs the Zp — 1 Gaussian lumi- 
nosity evolution curve at 15nm, while at 850nm there is a 
more clearer superposition of the 2 Gaussian scenarios. 




-3.5 -3 -2.5 -2 
Ig(Flux) {Jy} 



-1.5 -1 



Figure 6. Source counts of the high luminosity component only 
at 15|am (differential) and 850 |jm (integral) respectively. The 3 
model lines correspond to the three exponential density evolution 
scenarios of the ULIG component discussed in the text. Refer- 
ences for observed points are given in the text and figures. 



4 CONSTRUCTION OF NEW 
EVOLUTIONARY MODELS 

4.1 Source Counts 

From the analysis of the effects and constraints from density 
and luminosity evolution of the IR high luminosity popula- 
tion (assumed to be ULlGs) in the previous section, the 
following points / constraints can be made, 

(i) The evolution can be constrained by simultaneous con- 
sideration of the 15n-m and 850ti.m counts. 

(ii) The 15vim counts are predominantly dominated by 
low redshift galaxies and insensitive to the high redshift pop- 
ulation. 

(iii) The faint end of the SSO^im counts constrains the 
amount of density evolution allowed. 

(iv) Simple density evolution scenarios can fit the 15n.m 
counts but are unable to account for the sub-mm counts at 
fluxes > 2m Jy. 

(v) The addition of luminosity evolution has the effect 
of moving the source counts towards brighter fluxes in the 
lgiN)/lg{S) plane. 



expected, the addition of luminosity evolution moves the 
counts to the right in the lg{N) /lg{S) plane. 

At 15p.m, evolution at higher redshift {zp = 2) is not 
well constrained by the counts as the 15fi.m sources are ex- 
pected in general to reside at lower redshift. The Zp — 1 
Gaussian and plateau (due to the z~i contribution) scenar- 
ios, strongly violate the differential counts at IS^im not in 
fact due to the strength of the luminosity evolution alone 
but more due to the superposition of both the luminosity 
and density evolution effects at these redshifts. The 15^rm 
differential counts put a severe constraint on the magnitude 
of the luminosity evolution at lower redshift in the presence 
of density evolution. 

At 850nm, the effect of the addition of luminosity evo- 
lution is quite spectacular. The shift to the right in the 
lg{N) /lg{S) plane can effectively fit the SCUBA counts at 
fluxes both fainter and brighter than 2mJy. There is little 
difference between the high or low redshift evolutionary sce- 
narios implying that many of the SCUBA sources are at 



(vi) The 15 vim counts place strong constraints on the lu- 
minosity evolution at low redshift but are insensitive to lu- 
minosity evolution at higher redshift . 

(vii) Luminosity evolution at either high or low redshift 
can fit the 850[i.m counts (i.e. the SCUBA counts effectively 
only constrain the magnitude of the luminosity evolution in 
this scenario). 

From these considerations a model can be envisioned 
that includes both density evolution at lower redshift (peak- 
ing at 2 ^ 1) and luminosity evolution peaking at higher 
redshift. The exact details of the evolution are constrained 
by fits to the source counts which will depend on both the 
magnitude of the evolution and the degree of overlap in the 
redshift space between the tails of the density and luminos- 
ity evolution. The magnitude of the evolution at high red- 
shift will ultimately be constrained by the cosmic infrared 
background. 
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Figure 7. Source counts of the high luminosity component only at 15|am (differential) and 850^m (integral) respectively. The 3 model 
lines correspond to the three exponential luminosity evolution scenarios discussed in the text. Note that the low redshift density evolution 
scenario is also assumed and included in these plots. References for observed points are given in the text and figures. 



Using these 3 constraints (15vim differential counts, 
850 nm integral counts and the CIRB) and the results from 
the previous analysis, the following evolutionary scenario 
is presented. The basic framework of the cppRR model is 
retained with the galaxy population consisting of 4 main 
components. Normal (cirrus) galaxies are represented by 
the cool 60p.m IRAS luminosity function and are essentially 
non-evolving although a gradual decline towards higher red- 
shift is included to represent gradual formation and/or tran- 
sition through a starburst phase. The warm IRAS galax- 
ies (starbursts) and Seyfert galaxies are represented by the 
60|J.m and 12^tm IRAS luminosity functions respectively and 
evolve in luminosity at a rate oc (1 -|- z)^''^ to z=2.5 and then 
gradually decline at higher redshifts thus eliminating the 
artificial cut off employed previously in th e cppRR model 
and pointed out by Rowan- Robinson (200C) as being a gen- 
eral shortcoming with many models employing the backward 
evolution approach to galaxy evolution. Choosing a peak in 
the luminosity evolution at high redshift (z=2.5) is consis- 
tent with the higher redshift galaxies contrib uting more to 



the coun ts /background at longer wavelengths ( Gispert et al ^ 
200c|r7|The high luminosity end of the IR 60|i.m luminosity 
function, the LIG/ULIG component, lgL*sotJ.m ~ 11.61/q 



evolves in both number density and luminosity. The density 
evolution is strong and exponential, following eqn. ^ peak- 
ing at a redshift of 1 with the magnitude of the evolution, 
g — 220 and Gaussian width, a — 0.25. The luminosity evo- 
lution rises exponentially as eqn. ^ with a magnitude k = 40, 
a — 0.58 to a maximum redshift, Zp = 2.5 and then slowly 
declines exponentially to higher redshift. The magnitude of 
the luminosity evolution at the peak redshift in the star- 
burst, Seyfert and ULIG components is similar although the 
latter suffers a steeper decline towards lower redshift reflect- 
ing the more violent star formation within these galaxies. 

Fig. 1^ shows the fits of the new evolutionary model 
to the source counts as observed by ISO (15, 90, 170jim), 
IRAS (60^m) & SCUBA (850^m). The new model incorpo- 
rating the strongly evolving LIG/ULIG component fits the 
observed counts extremely well with the LIG/ULIG com- 
ponent accounting for the bumps in the 15nm counts at 
0.1-lmJy, the rise towards fainter fiuxes in the 60(i.m differ- 
ential counts, the excesses at 90(i.m & 170(i.m and much of 
the SCUBA population b righte r than 2mJy at 850(i.m. At 
ITOjim Matsuhara et al. (2000) concluded on the basis of 
the 170nm/90(i.m colours of galaxies in their Lockman Hole 
field, that the faint sources could not be cirrus dominated 
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Figure 8. Integral and differential counts at 15|-im, 60|am, 90|am 170|^m &; 850nm for the new evolutionary model of normal, starburst, 
ultraluminous and Seyfert galaxies as described in the text. Data is from IS^mi - Oliver et al. (1997), Aussel et al. (1999), Serjeant et al. 
(2000), Elbaz et al. (1999). 60^m - Hacking & Houck (1987), Rowan-Robinson et al. (1990), Saunders (1990), Gregorich et al. (1995). 
90iam - Matsuhara et al. (2000) (fluctuation analysis - dotted box) , Efstathiou et al. (2000). 170iam - Kawara et al. (1998), Puget et al. 
(1999), Matsuhara et al. (2000) (fluctuation analysis - dotted box). 850|am - Small et al. (1997), Hughes et al. (1998), Eales et al. (1999), 
Barger et al. (1999), Blain et al. (1999). 
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Figure 9. Number-redshift distributions at 15|am for all compo- 
nents in the model for fluxes ranging from 0.1-2mjy. The peak at 
the lowf st redshifts is due to the normal galaxy population. The 
ULIG p jpulation dominates at z~l. 



(normal galaxies) hut had to he starformirig systems, mnsis- 
tent with the model prediction of a change in the dominant 
population, from cirrus to star forming ULIGs at ~400mJy. 
Note that there is still significant discrepancy between the 
90[i.m ISO source coun ts co nducted by the the Japanese 
team (Matsuha ra et al. ( 200C ), Kawara et al. ( 200C )) and the 
ELAIS survey f|Efstathiou et al. 20nr ). The model predicts 



the 



a value lconservativcly between these 2 estimations. 



4.2 Number Redshift Distributions 

Figs. ^, 10, show number redshift distributions for the 
new evolutionary model at 15, 170 & 850ii-m respectively. 

The numerous ISO surveys at 15nm have a general con- 
sensus with a strong evolution of sources to 2; ~ 1. The cor- 
resp onding redshift dis tribution peaking in the range z~0.7- 
1.0 ( [Elbaz et al. 1999b| ). The N-z distributions in figure | are 
broadly consistent with this strong evolution although due 
to the form of the evolution assumed there is a double peak 
in the distribution. The lower redshift peak is due to the nor- 
mal galaxies that dominate the source counts to «lmjy. As 
the N-z distribution probes to fainter fluxes the contribution 
from the ULIG population (centred around z~l) increases. 
These models serve to show the effect of a strongly evolving 
population to z~l on the source counts and moving the peak 
of the density evolution to lower redshift - or increasing the 
Gaussian width of the evolution would produce a smoother 
N-z distribution effectively merging the 2 peaks. 

At 850(i.m, the redshift of the detected sources is notori- 
ously difficult to measure unambiguously, the main problem 
being the 2-3" (or more) positional uncertainty in t he sub- 
mm observations for a 15"resolution at 850nm (Hughes 
200c). Furthermore, there is often only a single definite de- 



tection at 850p.m in the entire spectral range from the sub- 
mm to the far-IR ~200(i.m . SCUBA preferentially selects 
galaxies at redshift > 1 due to the shape of the SED of 
galaxies in the sub-mm regime. A spectroscopic sur vey by 
Barger et al. (1999b) of a lensed sub-mm sample (smaiJ 



Figure 10. 

ponents in 
peak at low 
the N-z distitbution 
peak due to 
is narrow. 



Number-redshift distributions at 170|am for all com- 
model for fluxes ranging from 100-300mjy. The 
redshifts ~0.1 is due to normal galaxies. Note that 
is measured as dN/dz, therefore although the 
normal galaxies is high, the width of its distribution 



et al. 1999) concluded that the sub-mm population resided 
mainly between z«l~3 with a median redshift between 
1.5^^2. For the 3 SCUBA galaxies for which t here are CO 



spect r oscopy confirmed r e dshifts, all lie at z> l ( Frayer et al _ 
199^ ),( pVayer et al. 1999| ),( ^oucail et al. 199£| ). However, the 



first identifications from the Canada-U.K. deep sub-mm sur- 
vey have shown that the SEDs of the sub-mm sources (0.8, 
15, 450, 850vim and 6cm) are consistent with high IR lumi- 
nosity Arp 220 type luminous IR galaxies, many undergoing 
mergers/interactions with « 1/3 of the sample at 2: < 1, 
1 < z < 3 fc 2: > 3 an d a slightly high er median redshift 
of ~2 ( |Eales et al. 199c| ), ( |Lilly al. 1999| ). Eales et al ( |200C| ) 
combined th e Canada-UK deep survey, lensed surveys and 
HFF survey (Barger et al. 200C), producing a combined me- 



dian redshift of 2.41 
method of Carilli fc Yun 
Clements 



Finall y Sm all et al . (200C) using the 
([1999D, ([2OOOD (see also Dunne, 



Eales ( 2000 )) combined radio data for a sub- 
mm sample to fit radio-sub-mm SEDs thus obtaining rough 
estimates of redshifts finding a median redshift between 2-3 
for the sub-mm sources brighter than ImJy. In general the 
trend is for a spread of redshifts between 1^3 and higher 
with a median around z^2. Figure |ll| is broadly consistent 
with these estimates. At the brightest fluxes (~6mJy) the 
redshift-distribution is dominated by the ULIG galaxy pop- 
ulation consistent with the SCUBA observations of bright 
sources. In this region the evolution is dominated by the 
higher redshift luminosity evolution component. As the red- 
shift distribution is sampled towards lower fluxes two points 
are worth noting. The first is the emergence of the starburst 
galaxy population as the dominant contributor at high red- 
shift and the second is the strengthening of the number evo- 
lution component of the ULIG population. 
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Figure 11. Number-redshift distributions at 850 for individual components in the model for fluxes ranging from l-6mjy. In general, 
the normal galaxies are represented by the low redshift peak in the N-z distribution. 



4.3 Background Radiation 



expected to com e from relatively low redshift galaxies (Gis- 



Further constraints on the evolutionary model come from 
the cosmic infrared background (CIRB) which manifests it- 
self as the superposition of the unresolved population of 
faint sources in the Universe. The far-infrared background 
has been well constrained from 100/im — 2mm by mea- 
s urements using th e F IRAS /DIRBE instru m ents on CQBE 
dPuget et al 1996| ), ( [Fixsen et al. 199^ ), jHauser et~ar 
199q)a id from ISO sourc e c ount fluctuation analysis at 240- 



140^ ^m (iPuget et al 1999| ), ( |Matsuhara 2000| ). Finkbeiner et 
al. ( 200(]| ) have also reported possible detections and excess 
at 60 & lOOvim. In the sub-mm, lensed counts to 0.5mJy 
have discovered enough sources to approximately account 
f or the SSOyim bac kground of !/J^ = 5 x 10~^nWm~^sr~^ 
(Blain et al. 1999a). In fig. |l2| the contribution of the new 
evolutionary model to the CIRB is shown compared with 
the observations described above. 

The model lies comfortably within the upper and lower 
constraints set by FIRAS. The cosmic backgroun d pea ks at 
~140ii-m where detection s by both Hauser et a l. (1998) and 
the FIRBACK survey ( Lagache et al. 199E| ) are in good 
agreement. Note that the longer wavelength (240^.m) mea- 
surements of COBE a nd FI RBACK do not agree so well al- 
though Lagache et al (1999) assumed a significant correction 



due to warm interstellar dust (WIM) in their calculations. 
The model matches both the position and the magnitude 
of this peak extremely accurately. One direct consequence 
of this is that the majority of this peak emission would be 



pert et al. 2000). The evolution in the high redshift regime is 



strongly constrained by the CIRB/Sub-mm background at 
longer wavelengths where the CIRB slope is shallower than 
the average galaxy SED implying a strong contribution from 
the integrated light of high redshift sources where the peak 
of the IR-SED is being sampled in the sub-mm waveband. 
It should be noted that assuming an O = 1 cosmology re- 
sults in an extremely mariginal fit to the longer wavelength 
(sub-mm - mm) regime. 

Constraints on t he CIRB at shorter wa velengths come 
from DIRBE, IRTS ( iMatsumoto et al. 200c|) , galaxy source 
counts and TeV gamma rays ( Stecker 200(]| ), although the 
source count measurements and background measurements 
are yet to converge on an agreed value probably due to the 
difflculty in ascertaining the correction fact or due to zodiacal 
light that peaks at ~25^m (Dzernoy 2QQC). 

Although the peak of the CIRB is well explained by 
the strongly evolving ULIG population in this model, the 
stronger contributor at sub-mm wavelengths are in fact the 
starburst galaxies. Furthermore although much of the sub- 
mm background may be made up from a population of ULIG 
galaxies (30-50% being resolved into individual galaxies with 
Ss50 > "imjy), the surfac e density of t he sub-mm popula- 
tion is_significantly higher (Hughes 2000). Note that Eales et 



al. (2000) have pointed out that the SCUBA fluxes may be 
upwardly biased by a factor of 1.4 due to source confusion 
and noise, resulting in a reduction of the resolved CIRB at 
850^m to 20%. 
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Figure 12. The contributions to the infrared-sub-mm background light predicted by the model. The Shaded area in the figure represents 
the constraints set by the FIRAS experiment on COBE. Source Counts - refers to lower limits set from the integrated light of galaxy 
source counts. Gamma Ray - refers to limits set by the attenuation of TeV gamma rays by the background radiation. FIRBACK - the 
integrated light of galaxies calculated from the FIRBACK survey (Puget et al. (1999)). 



Such models of the CIRB are also consistent with the 
850ii-m source counts (see fig. |^. The ULIG component is 
completely dominant at bright fluxes down to ~6mJy. Be- 
low this limit the starburst galaxies begin to contribute more 
strongly, where at Sgso < 2mJy the y bec ome the dominant 
source population. Peacock et al. ( ^00(]| ) made a statisti- 
cal analysis of the 850(i.m SCUBA map of the Hubble Deep 
Field (HDF) estimating that as much as 48% of the CIRB 
at 850 ^tm could come from the sub-mm emission f rom the 



UV st arfor ming galaxies with 30% from U LIG's (Hughes 



et al. 



199S) 



J398|) and a further 10% from AGN ( [Almaini et al 
f the starburst population in the model is literally 
taken as the equivalent of the UV galaxies in the HDF then 
the model predictions slightly over predict this contribution 
from the starburst galaxies although the transition between 
starburst-LIG-ULIG is not clearly defined and therefore has 
a large scope for interpretation. The sub-mm source count s 
are also consistent with the results of Chapman et al. (200C) 
who carried out a sub-m m survey of 33 Lyman break galax- 
ies (Madau et al. 1996) to ~1.3mJy. The average fiux of 
these sources was found to be Sgso ~ O.GmJy implying that 
they would occupy the fainter end of the sub-mm source 
counts and would not significantly contribute to the brighter 
Sgso ~ 5mJy source counts. Therefore it would seem that at 



the brighter sub-mm fiuxes (2-5mJy) we are observing the 
ULIG population where as much as one third of the CIRB at 
850 itm can be resolved into individual bright sources. The 
underlying fainter population are distinct from the ULIGs 
and are more akin to the UV Lyman Break Galaxies and / or 
IRAS starburst galaxies. In this scenario the 850itm counts 
are a superposition of 2 galaxy populations, starbursts and 
ULIGs. 



5 DISCUSSION 

The aim of this work has been to find an evolutionary model 
that can viably fit the galaxy source counts from the mid- 
IR to sub-mm wavelengths whilst not violating the con- 
straints set by the IR background and CMB measurements 
and was not born from a fundemental physical assumption. 
The predicament of the model discussed in this paper is how 
does the bimodal evolution (i.e. dual evolution) manifest it- 
self in the evolutionary history of a galaxy. Are we seeing a 
representation of two stages of galaxy evolution of the same 
population or different evolutionary paths of separate pop- 
ulations? Of the 19 SCUBA sources detected at 850iJ.m in 
the Canada-UK sub-mm survey, only 2 were found to have 
ISO 15nm associations implying that the ISO and SCUBA 
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sources are in fact different po pulations (or different red- 
shift regimes) (Eales et al. 2000). Furthermore, what role do 



the ULlGs play in the evolutionary history of AGN? Do all 
QSOs go through a dusty ULIG phase seeing that the space 
density and luminosities of QSOs and ULlG s in the local 
Universe are quite similar (Soifer et al. 1987). 



If this scenario is considered as 2 phases of evolution 
then the initial formation of the core (collapse of dark mat- 
ter halo) would take place at high redshift accompanied by 
the initial formation of a black hole and associated star for- 
mation slowly declining towards lower redshift as the initial 
fuel is used up. The peak in the co-moving number density 
of AGN is at redshifts ~2, at lower redshifts these AGN 
are fou nd to be hosted in giant elliptical galaxies ( Mc- 
Clure e |t al. 199c| ). Thus the mid-IR 3-30um SED may be 



indicative of the galaxies merging at lower redshifts acti- 
vating / re-activating the central black hole where at some 
time the buried AGN expels the surrounding dust becom- 
ing the dominant energy source. The SED >50p.m remains 
starburst dominated. The ULIGs would undergo rapid evo- 
lution from high redshift eventually evolving into today's 
local giant elliptical galaxies. Local giant ellipticals (gE 
galaxies) are the only local systems with core densities high 
enough {100Mqpc~^) to account for a pa st ULIG evolution- 
ary phase (Korm endy & Sanders (1992), Binggeli (1994), 
Trentham (|200C|)). 



6 CONCLUSIONS 

The backward evolution approach to galaxy evolution has 
been re-visited in light of the recent observations by ISO 
and SCUBA in the infra-red and sub-mm regions. Assum- 
ing classical evolutionary scenarios incorporating power law 
luminosity and density evolution it is impossible to simul- 
taneously fit all the source counts from 15-850iJ.m and the 
CIRB. 

A new evolutionary model is investigated which as- 
sumes the foundations of the Pearson & Rowan- Robinson 



(1996) model (cppRR). In addition to this original frame- 
work, a strongly evolving high luminosity IR galaxy com- 
ponent (referred to as the ULIG component but specifically 
the high luminosity end of the 60n-m luminosity function) 
is added. This new population evolves in both density and 
luminosity and is capable of fitting both the source counts 
from 15-850|i.m and the cosmic far-infrared background ra- 
diation, including the peak at ~140vim and the long wave- 
length tail. 

In this scenario the ULIG component is tenuously in- 
terpreted as the progenitor of the present day giant elliptical 
galaxy population. 

The excesses in the 15n.m and 170^i.m counts are well 
explained by strong evolution to redshift ~1 of a popula- 
tion of luminous, IgLia ^ 12Lq, infrared galaxies. This 
population of galaxies becomes prominent at fiux levels just 
below the limits of the IRAS 60vim survey (although previ- 
ously hinted at) and can be associated with the archetypical 
ULIG Arp220. These results are generally consiste nt with 
the numerous surveys at 15^.m ( Elbaz et al. 200C ) which 
found that the dominant population above a redshift of ~0.5 
were massive {M > 10"Mq, luminous Ljr > 10""^^I/ 
galaxies. 







At 850ii.m, these galaxies dominate the bright SCUBA 
counts at S^6mJy. Towards fainter fiuxes the starburst 
galaxies begin to contribute strongly to the source counts, 
becoming dominant at fiuxes < 2mJy. In this scenario 
the sub-mm counts can be divided into 2 distinct popula- 
tions, the dusty, giant ULIGs which would be considered 
the progenitors of giant elliptical galaxies and the IRAS 
starburst population / Lyman break galaxies. This conclu- 
sion is il l goo d agreement with the analysis of Busswell & 
Shanks (200C) who transformed the optical counts of faint 
starburst gal axies to the sub-mm waveb and. Using a Bruzual 
& Chariot (Bruzual & Chariot 1993) galaxy model with 
an exponentially decaying SFR of r = 9Gyr, their models 
successfully pre dicted optical result s, including the Lyman 
break galaxies (Metcalfe et al. 200C). Furthermore they ac- 
counted for the sub-mm counts below a sensitivity of ~lmjy 
but failed to fit the counts at brighter fiuxes, concluding that 
an additional population was required (ULIGs/AGN?). Of 
further interest is that their models were also deficient in 
fitting the CIRB at 140-240^m where the ULIG population 
is most prominent. 

Note, analysis of the so called Madau plot of the star 
formation history of the galaxies in this evolutionary model 
will be addressed in a later paper. The primary aim of this 
work is to ascertain what evolutionary scenarios can succes- 
fuUy fit the observed source counts / background radiation. 

From an analysis of luminosity and density evolution 
the following conclusions can be drawn from the galaxy 
source counts. 

(i) The mid-IR region (15ii.m) is sensitive to evolution in 
the low redshift population but is not strongly constrained 
by evolution at higher redshift (z>l) 

(ii) The far-IR region at 60-90(i.m does not presently in- 
clude enough information to constrain the evolution in the 
source counts. 

(iii) The far-IR region at 170^.m is sensitive to the low 
redshift population and the faint end counts constrain the 
contribution to the peak of the CIRB. 

(iv) The sub-mm counts constrain the magnitude of the 
evolution in the high redshift Universe z>l but due to the 
steep K-corrections, cannot distinguish evolutionary forms 
between redshifts 1-2. 

Using the short wavelength 15(i.m and long wavelength 
850ii-m counts, together with the CIRB as bounding condi- 
tions, the evolution over a wide range of both wavelength 
and redshift can be constrained. It is found that a scenario 
where the ULIG population undergoes 2 phases of evolution, 
a rapid merging epoch to z~l and exponential evolution in 
luminosity to higher redshifts, best fits the observations. 

In this scenario the high redshift ULIGs detected at 
850 vim fiuxes J5 2m Jy may be interpreted as the progenitors 
of todays giant elliptical galaxies while the population at 
fiuxes fainter than Ssso ^2mJy represent the IRAS starburst 
population / Lyman break galaxies. 

Obs ervations in the next few years in the mid-IR with 
SIRTF ( [Rieke 2000]) an d the Japanese ASTR O-F mission 
(M atsuha ra et al. ( |l99l ), Watarai et al. ( |2000| ), Pearsonet 
al. (I2OOII)) , in the far-IR-sub-mm -mm with FIRST ( [Pilbrat 



G.L. 2OOOD , LMT dSchloerb 1997| ) and an all sky sur vey a t 



far-IR wavelengths with ASTRO-F (Kawada et al. ( [1998D 
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Takahashi et al. (2000)) will undoubtebly solve many of 
these mysteries. 



7 INTERNET ACCESS TO COUNT MODELS 

The source count models can be accessed through the world 
wide web at 



htt : / /www.ir.'i 



.jp/ ^ cpp/ counts/ 
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